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The LEP experiments give a lower bound on the neutrahno mass of about 46 GeV which, however, 
rehes on a supersymmetric grand unification relation. Dropping this assumption, the experimental 
lower bound on the neutralino mass vanishes completely. Recent analyses suggest, however, that in 
the minimal supersymmetric standard model (MSSM), a light neutralino dark matter candidate has 
a lower bound on its mass of about 7 GeV. In light of this, we investigate the mass sensitivity at 
the ILC for very light neutralinos. We study slepton pair production, followed by the decay of the 
sleptons to a lepton and the lightest neutralino. We find that the mass measurement accuracy for 
a few-GeV neutralino is around 2 GeV, or even less if the relevant slepton is sufficiently light. We 
thus conclude that the ILC can help verify or falsify the MSSM neutralino cold dark matter model 
even for very light neutralinos. 



I. INTRODUCTION 



The supersymmetric Standard Model (SSM) [H, is a 
well motivated extension of the Standard Model of par- 
ticle physics which solves the hierarchyproblem between 
the weak scale and the Planck scale [3 Q- In order to 
guarantee a stable proton, usually a discrete symmetry 
beyond the SM gauge symmetries is imposed which pro- 
hibits the baryon- and lepton-number violating terms in 
the superpotential: R-parity proton hexality Q, or 
a Z4 R-symmetry Q. This is then called the minimal 
supersymmetric Standard Model (MSSM). The discrete 
symmetry furthermore guarantees that the lightest su- 
persymmetric particle (LSP) is stable and thus a dark 
matter candidate If it is to constitute the en- 

tire dark matter in the universe it must be electrically 
and color neutral Q . Here we focus on the lightest neu- 
tralino Xi the LSP. In order to avoid overclosure of 
the universe, the neutralino must be either very light 
M^o < O(leV) (Cowsik-McLelland bound) [H, [H or 
heavy M^o > 0(10 GeV) (Lee- Weinberg bound) 
We shall make this latter lower bound more precise 
[2]| . In this paper we are interested in how a neutralino 
mass close to the Lee-Weinberg bound could be mea- 
sured at the ILC. This is potentially a stringent test of a 
MSSM light dark matter model. 
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Within the MSSM the spin-1/2 superpartners of the 
hypercharge B boson, the neutral SU(2) W boson and 
the two CP-even neutral Higgs bosons mix after elec- 
troweak symmetry breaking. The resulting four mass 
eigenstates are the neutralinos and are denoted x°, i = 
1, . . . , 4. The masses are ordered M^o < . . . < M^o. If pro- 
duced at colliders, the lightest neutralino behaves like a 
heavy stable neutrino and escapes detection. The spin- 
1/2 superpartners of the charged SU(2) W boson and of 
the charged Higgs boson also mix after electroweak sym- 
metry breaking. The resulting mass eigenstates are the 
charginos and denoted XiLi 21 with ordered masses. See 
Appendix A for details. 

The current Particle Data Group (PPG) mass bound 
from LEP on the lightest neutralino is |22h25| 

M^o > 46 GeV. (1) 

This bound is obtained by searching for charginos and 
thus setting a boimd on the SU(2) gaugino mass M2 and 
the Higgs mixing parameter /i. Using the supersymmet- 
ric grand unified theory (SUSY GUT) relation between 
M2 and the U(l)y gaugino mass term Mi 

Ml = ^ tan^ M2 , (2) 

the chargino search can be translated into a bound on 
Ml. The neutralino mass matrix is computed for all al- 
lowed values of the supersymmetric parameters, taking 
into account Eq. as well as the lower bound on the 
ratio of the Higgs vacuum expectation values, tan /3 > 2, 
from the LEP Higgs search [2^. Then one obtains as 
the lowest possible neutralino mass the bound in Eq. ([Ij . 
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If, however, the assumption Eq. ^ is dropped, there is 
no lower laboratory or astrophysical bound on the neu- 
tralino mass [l^ [26l - l32| . Even a massless neutralino is 
allowed. This is now included in the PDG as a comment 
[33| . There is however a cosmological bound which we 
now discuss. 

Relaxing the SUSY GUT assumption in Eq. ©, it 
is possible to derive the Lee-Weinberg lower limit on 
the mass of the neutrahno LSP, MT", in the MSSM 
with real parameters. It was first determined for large 
pseudoscalar Higgs masses [l^ [111, obtaining Af""" — 
0(15GcV). It was subsequently realized however [TtI . 
Il8| . that a region of parameter space exists with a low 
pseudoscalar Higgs mass and high tan /3, in which the 
neutralino lower mass limit reaches M""" w 6 GeV. This 

Xi ^ 

is due to an enhancement in the neutralino annihila- 
tion cross section from annihilation to b-quarks via Higgs 
bosons, which keeps the predicted relic density below the 
observed limits. This was confirmed in Ref. |34l |. There it 
was furthermore noted that this area of parameter space 
would be testable at the Tevatron, for example, with the 
Higgs search results in which the Higgs is produced in as- 
sociation with a b-quark, as well as via the Bf, — » tJ. '^ii'' 
limit. See also the more recent work in Ref. [35|-|38|. In 
a very recent paper [111, the authors argue that these 
constraints have a relatively minor impact on the light 
neutralino parameter space of the MSSM [s^, and that 
the lower bound is 



evidence for non- minimal supersymmctric models and/or 
non-standard cosmologies. 

The paper is organized as follows. In section|lIl we pro- 
vide an overview of some of the methods that have been 
suggested to measure the lightest neutralino mass at the 
ILC. We then focus, in Section HID on the neutralino 
mass measurement that can be done using slepton pair 
production, and make a first estimate of the accuracy of 
this method for a very light neutralino. In Section llVi 
we describe a simulation of slepton pair production at the 
ILC, and use it to make a better determination of the ac- 
curacy attainable for the neutralino mass measurement. 
We summarize and conclude in Section IVl 



II. NEUTRALINO MASS MEASUREMENTS 

Several methods have been suggested in the literature 
to measure the mass of the lightest neutralino at the ILC 
[64l - |68| . Throughout, the authors have focused on a neu- 
tralino heavier than the LEP bound in Eq. ([T]) . For exam- 
ple, the widely studied SPSla point (without a slope) [g^ 
has Af^o = 97.1 GcV. The most straightforward method 
involves considering slepton pair production, followed by 
the decay of each slepton to the lightest neutralino and 
a charged lepton 



(4) 



M"r w 7 - 8 GeV . 

Xi 



(3) 



Recently there has also been an increased interest in 
light dark matter candidates with a mass of order 5 GeV 
due to the DAMA/LIBRA and CoGeNT [U direct 
search results. Ref. [i^l suggests that the required scat- 
tering cross sections in the detectors cannot be obtained 
within the MSSM, though in Ref. ^ it is claimed that 
the current constraints do in fact allow sufficiently high 
cross sections for these experimental hints to be explained 
by an MSSM neutralino. In the NMSSM, it is similarly 
claimed in Ref. [H, that high cross sections cannot 
be obtained, though Refs. [H |4Q| identify regions of pa- 
rameter space in this model in which the presence of a 
light singlet Higgs can lead to large enough cross sec- 
tions. The authors of Refs. [131 argue that solutions also 
exist for an extended NMSSM. The experimental results 
in Refs. [H, E3 have lead to a flourish of alternative 
schemes, e.g. |48l - [58| . We note, however, that these mod- 
els are severely constrained by the CDMS- II fsol [60} and 
XENONIO ^ and XENONIOO data [6l[6|^ 

Given these considerations it is thus of great interest 
how well the neutralino mass can be determined in the 
low-mass region. It is the goal of this paper to test the 
sensitivity of the mass measurement at the ILC. There is 
the possibility that a neutralino mass will be measured 
which is too small to be reconciled with the observed 
relic abundance, if the real MSSM and standard cosmol- 
ogy arc assumed. Such a measurement would be striking 



Here the (s)lcptons are restricted to the first two genera- 
tions. The measurement via the third generation (s)tau is 
diluted by the additional decay to the neutrino(s). Mea- 
suring the energies of the final state leptons, one can ex- 
tract information on the neutralino and slepton masses. 
The typical relative precision achieved is in the per mille 
range [13, l66l - [68| . We go beyond this work and discuss 
this method in detail for a very light neutralino. 

A second method in the literature is based on the pair 
production of the second lightest neutralinos. This is 
followed by the decay of each neutralino via a (virtual) 
slepton to a charged lepton pair and the lightest neu- 
tralino. 



A2A2 



(5) 



where each X2 decays independently and thus £1 need 
not equal £2- In fact, the case £1 # £2 reduces the com- 
binatorial uncertainty. In Refs. [TO, UM ^he authors then 
propose to measure the di-lepton invariant mass and the 
di-lcpton energy and to use these to measure the two 
lightest neutralino masses. 

Of necessity these methods also always involve other 
supersymmctric particles and their masses. For exam- 
ple, the first method relies on the production of sleptons. 
The second method relies on the production of the second 
lightest neutralino and then its decay to an intermediate 
slepton. Thus both of these methods can be improved 
by measuring the corresponding supersymmctric masses 
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directly. For example the slepton mass can be well deter- 
mined by an energy scan over the slepton mass thresh- 
old [7^. Similarly a scan over the production threshold 
energy of the process given in Eq. ((S]) gives a tight con- 
straint on the mass M^o [t^] ■ 



III. SLEPTON PAIR PRODUCTION AND THE 
NEUTRALINO MASS 

In this section, we study the measurement of the light- 
est neutralino mass using slepton pair production at the 
ILC, as shown in Eq. (HJ. The slepton decay to a lepton 
and the lightest neutralino is a two-body decay. There- 
fore in the slepton rest-frame the lepton energy is com- 
pletely fixed by the slepton and neutralino mass. Ig- 
noring initial and final state radiation (ISR and FSR), 
beamstrahlung, and detector effects for the moment, the 
slepton energy is then just the beam energy. Thus the 
lepton's lab-frame energy Ei is fully determined by the 
angle with which the slepton emits the lepton in the 
slepton rest-frame. The angle is measured with respect 
to the slepton lab momentum direction. We then have 
for the lepton energy 



E,, 



A/? 



(l + /3cos6io) 



(6) 



Here /3 = y 1 — 4m? /s is the slepton velocity in the lab 

frame, \/s/2 is the beam energy, and denotes the 
slepton mass. The event distribution of Ei is flat between 
its maximum Ej^, when cos^o = li a-nd its minimum 



E^, when cos6'o 



-1. The equations for E^ and E^ 



can be inverted to find the slepton and neutralino masses 
squared in terms of these endpoints. 



Mj = s 



E+E^ 



iE+ + E^y 



and 



Xi 



1 



E. 



E^ 



Vi/2 



(7) 



(8) 



We have listed the squared formulas for later use. Taking 
the positive square root we then obtain for the masses 



E,+E. 



and 



E, 



E. 



(9) 



(10) 



The sensitivity of the neutralino mass measurement thus 
depends on the accuracy with which E± can be measured. 
Looking at Eq. ©, it is clear that E± only have a 



weak dependence on AI^ 



for M^o < M^. 



Thus we would 



expect the accuracy of the neutralino mass measurement 
to deteriorate for sufficiently small neutralino masses. We 
set out to quantify this below. 

Limited statistics and detector and beam effects in- 
troduce uncertainty into the endpoint determination. 
Nonetheless, for typical slepton and heavy neutralino 
masses, the endpoints and the masses can be determined 
to sub-GeV accuracy For very light neutralinos, 

however, even small errors in the endpoint measurements 
can lead to a large fractional error in the neutralino mass 
determination. 

Before studying this issue with a simulation, we can es- 
timate the mass determination accuracy for a light neu- 
tralino by combining the quoted accuracy from an exper- 
imental study by Martyn [gtI ] with a simple error analy- 
sis. Assuming that E-^- and i?_ are independent random 
variables, then from Ec^. pOI) we can derive 



M? 



SE^ 



SE^ 



(11) 



For light neutralinos, the first term in Eq. (jll[) is negligi- 
ble. The second term dominates and is identical for E^ 
and £'_, so we can write 



SAEo 



M? 



5El 



SE^ 



(12) 



In the simulation we describe below, we consider SUSY 
scenarios with Mg^ = 100 and 200 GeV and varying 
M^o, a center-of-mass energy y/s = 500 GeV and an 
integrated luminosity £ = 250 fb^^. For illustration, we 
here assume these experimental parameters as well as a 2 
GeV neutralino mass and 100 GeV selectron mass. Thus 
the factor in front of the square root in Eq. (fT2)) is 10. 

In Ref. [13 , the error on the endpoint determinations is 
given as 5E+ = 0.11 GeV and SE_ = 0.02 GeV for a sce- 
nario with M^o 



93 GeV, Mi = 143 GeV, C = 200 fb" 



and y/s = 400 GeV. In this scenario, because the neu- 
tralino is heavy the two terms in Eq. (|lip are comparable, 
so we cannot use Eq. (|12p . Using Eq. ([TlT) instead, we 
obtain SPEo ~ 100 MeV, which agrees exactly with the 
quoted result of the detailed study in Ref. J67[. 

To translate these into an estimate for SE± in our sce- 
nario, we need to take into account several modifications, 
(i) The endpoint locations have changed significantly be- 
cause the slepton and neutralino masses are different. 
Therefore also the experimental energy resolution at the 
endpoint locations is different, (ii) The number of events 
for slepton pair production is different in our scenario 
due to the different masses, center-of-mass energy, and 
luminosity, so that the statistical error on the endpoint 
determination is different. Determining the effect of (i) 
requires choosing a parametrization of the detector's en- 
ergy resolution, which we discuss in the next section and 
provide in Eqs. (fT4)) and ([T5|). Taking these two factors 
into account, we can estimate the ratio of our endpoint 
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energy uncertainty to Martyn's in Ref. [67 



SE]i 



Martyn 



j^Martyn 



SEcxp{E = E. 



Martyn 
± 



(13) 



where we have estimated that the uncertainty in the end- 
point determination drops with the square root of the 
number of observed events. Plugging the relevant num- 
bers into the above expression for which dominates 

the error, in fact yields ^^Mattyn — 1-3, since the increased 

number of events in our scenario is partially canceled by 
the reduced detector resolution at the higher value of E^^ . 

Referring again to Eq. (jl2p , we can then estimate that 
SM^o ~ 1.4 GeV for our scenario. In other words, for a 
neutralino mass of 2 GcV, the mass can be determined 
to about 70% accuracy. This suggests that a useful mass 
measurement can be performed for very light neutralinos, 
and in particular in the range M^o ~ 5 GeV that is 
particularly interesting for dark matter phenomenology, 
sub-GeV accuracy should be possible. 

On the other hand, if we carry out the same estimate 
for a 2 GcV neutralino and instead a 100 GeV en, we find 
that the factor in front of the square root in Eq. is 
now 40, and the cross section for selectron pair produc- 
tion is also lower so that the statistical uncertainty is 
larger. In this case we find SM^o ~ 15 GeV, suggesting 
that in this case at best an upper limit on the neutralino 
mass can be set. 

While this simple estimate gives a qualitative illustra- 
tion of the difficulty of measuring a light neutralino mass, 
we would would like to check it with a more thorough 
analysis and more precisely quantify the accuracy possi- 
ble for a light neutralino mass measurement at the ILC. 
We do this in the next section. 



IV. SIMULATION OF NEUTRALINO MASS 
MEASUREMENT FROM SLEPTON PAIR 
PRODUCTION 

Thanks to the simple kinematics of slepton pair pro- 
duction, it is possible to estimate the precision for a 
Xi mass measurement at the ILC from a rather simple 
Monte Carlo simulation. We describe this in the follow- 
ing. 

First, the number of produced slepton pairs for a given 
centre-of-mass energy ^/s and luminosity C is calculated 
for a beam polarisation of {'P^,'Pe+) = (+80%, -60%) 
using the program SPheno [72| . which implements the 
cross section formulae from Refs. [tsI - Itg} . This choice of 
signs for the beam polarization maximizes the produc- 
tion cross-section. For each event, two lepton energies 
are thrown according to a flat probability density dis- 
tribution between E- and E+. In order to take effects 
caused by beamstrahlung into account, E- and E^ are 
evaluated for each event using the reduced centre-of-mass 



spectrum computed by GUINEA PIG [77[. The energy 
difference ^/s — ^/s' is lost in the form of beamstrahlung 
photons. As a result the sharp edges in the lepton energy 
spectrum are smoothed out a bit. 

The resulting lepton energies are subsequently smeared 
according to the expected momentum and energy resolu- 
tion. This smoothes out the edges even further. For elec- 
trons, the minimum of track momentum resolution and 
the energy resolution of the electromagnetic calorimeter 
(ECAL) for the considered electron energy is employed. 
In the case of muons, the momentum resolution of the 
tracking system is always used. For these quantities, the 
following parametrizations are used: 



(14) 



A— = 1 • lO^"^ GcV"^ (tracker), 
Pt 

.0.011 (ECAL). (15) 



E 



vAB/GeV 



Any polar angle dependence of the tracker resolution is 
neglected. Instead a rather conservative average reso- 
lution is applied (compare e.g. Ref. [ill). We checked 
that the results do not depend strongly on the assumed 
tracker resolution since for the considered SUSY masses, 
the Xi mass measurement is dominated by the calorime- 
ter resolution. The above parametrization of the ECAL 
resolution which wc employ is the one obtained with a 
detector prototype in test beam measurements [79| . The 
effects of a limited detector acceptance, signal selection 
cuts and inefficiencies in the electron and muon recon- 
struction are approximately accounted for by applying 
an overall efficiency of 50%. This roughly corresponds to 
the values obtained in Ref. ^6T\ using a more detailed sim- 
ulation. This more detailed study also showed that back- 
ground rates are rather small [70[. Therefore, outside of 
this overall efficiency, we neglect backgrounds completely 
in our study. 

The edge positions of the lepton spectrum obtained in 
the described way are finally fitted using an unbinned 
likelihood fit. The fitted shapes are 




E <E^ 
E> E^ 



(16) 



UiE) 



^erfc 
ierfc 



( V2f+ ) 



\/2cr 



E <E. 



+ 



E>E+ 



(17) 



for E^. If one chooses af =0"^, Eqs and pT|) are 
the results of a convolution of an upward and a down- 
ward step function with a Gaussian. Between the nom- 
inal edge positions E^ and E^ the shape of the lepton 
energy spectrum is influenced by beamstrahlung and en- 
ergy/momentum resolution, whereas outside the nominal 
edge positions, the shape is only determined by the en- 



energy which is thrown according to the luminosity ergy/momentum resolution. For this reason af^ and cr- 
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FIG. 1: Estimated precision of the Xi raass measurement 
from e-RCR production as function of the Xi rnass for en 
masses of 100 GeV and 200 GeV. The yellow bands repre- 
sent the estimated uncertainty of 30% related to the sim- 
plifications of the Monte Carlo simulation used. The as- 
sumed centre-of-mass energy is ^/s = 500 GeV, the inte- 
grated luminosity C = 250 fb~^ and the beam polarization 

{'Pe-,'Pe+) = ( + 80%, -60%). 



are treated as separate parameters in the fit. The fitted 
values of the parameters E- and i?+ do not in general 
coincide with the values of E- and E+. The reason is 
that the asymmetric shape of the beamstrahhmg energy 
spectrum leads to a certain offset. To correct for this 
bias, a Monte Carlo based calibration procedure is used. 

The uncertainty on the edge positions, and thus the 
masses, is determined by creating toy Monte Carlo 
datasets. For each toy data set, the fitted and corrected 
values of E^ and E- can be converted into the squared 



neutralino mass using Eq. 



The distribution of 



from the ensemble toy data sets is approximately Gaus- 
sian, and we use the width of this distribution to deter- 
mine the uncertainty on the mass measurement. For low 
neutralino masses, the distribution can have support in 
the unphysical region where rn^o < 0. To account for 

this, we used the Feldman-Cousins method [s^l to have 
a smooth transition between a mass measurement, which 
is possible for heavier neutralinos, and an upper bound, 
which is necessary for very light neutralinos. 

Using the described procedure, we agree within 
roughly 30% with the results in Refs. [6^[63,[23], which 
were obtained for M^o = 71.9, 96, and 135 GeV, respec- 
tively. Therefore we assign a systematic uncertainty of 
30% to our results due to the simplifications of our sim- 
ulation. 

Our estimate of the precision of the Xi mass mea- 
surement from ejieji production at the ILC is shown 
in Fig. [T] as a function of the Xi mass. The assumed 
luminosity is 250 fb~^ at a centre-of-mass energy of 
y/s = 500 GeV with a beam polarization of (Ve- , 'Pe+ ) = 
(+80%, —60%). Even for Xi masses as small as 2 GeV, a 
precision on the Xi mass measurement of ~ 0.6 GeV can 



We find that below about 2 (4) GeV for a 100 
(200) GeV selectron, a mass measurement is no longer 
possible and we can only set an upper bound on the neu- 
tralino mass. For example, for m^o = 1 GeV, the 95 % 
CL upper limits are 2.5 GeV (7.6 GeV) for a selectron 
mass of 100 GcV (200 GeV). 

We note that the precision of the mass measurement 
that we obtain in this simulation is roughly a factor of 
two better than the rough estimate of the precision made 
in Section IIIII This is due to the simplistic scaling as- 
sumptions made there, namely that the endpoint energy 
determination accuracies scale like as the number 

of events changes, and linearly with the detector reso- 
lution as the endpoint energy changes. Using dedicated 
simulations we find some deviation from this simple scal- 
ing which is due in a large part to the effect of beam- 
strahlung. This more realistic scaling can account for 
the discrepancy between our estimate and simulation re- 
sults. 

Combining the results from PlrJir production with 
those from eneu production does not lead to a sizable 
improvement of the obtained precision. The reason is 
the significantly higher cross-section for e^e^ produc- 
tion due to the additional i-channel contribution, which 
is especially important for low neutralino masses, lead- 
ing to a factor of 2 to 3 weaker constraints from jinP-R 
production. 



V. SUMMARY AND CONCLUSION 

A light neutralino in the several- GeV mass range is cur- 
rently of special phenomenological interest. Recent dark 
matter direct detection experiments hint at the possible 
existence of such a light particle. On the other hand, 
recent phenomenological analyses claim that an MSSM 
light neutralino dark matter candidate has a lower bound 
on its mass around 7 GeV. 

If a light neutralino exists, it would therefore be ex- 
tremely important to obtain an accurate determination 
of its mass. Techniques for measuring neutralino masses 
at the ILC have been developed and shown to have ex- 
traordinary precision for the more conventional 50-100 
GcV range. These techniques, however, have not been 
studied for much lighter neutralinos. 

In this paper, we have studied one of these 
techniques — measuring the lepton energy spectrum in 
slepton pair production events — and determined its use- 
fulness for the measurement of very light neutralino 
masses. We showed with a simulation that this technique 
continues to have useful accuracy for a neutralino with 
a mass as low as a few GeV. For example, we showed 
that it is possible to measure the mass of even a 2 GeV 
neutralino to sub-GeV accuracy if the mass of the right- 
handed selectron is 100 GeV. For a 200 GeV selectron, 
the precision is about 2.7 GeV for a 4 GeV neutralino. 
For even lighter neutralinos, we showed that this method 
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can give an 95 % CL upper bound of 2.5 (7.6) GeV for a 
100 (200) GeV selectron. 

Such mass measurements at the ILC will thus be indis- 
pensable in testing the MSSM thermal cold dark matter 
picture if a very light neutralino exists. 



Appendix A: Chargino and Neutralino Mixing 

Here we summarize the mixing of the electrowcak gaug- 
inos and Higgsinos, which we use in the paper. The spin- 
1 /2 superpartners of the gauge bosons and the scalar 
charged Higgs field, mix after electroweak symmetry 
breaking. The resultin g rn ixing matrix in the wino, Hig- 
gsino basis is given by [8l[ 



il/2 V^Mw sin /3 
\/2Mn/ cos/3 fj, 



(Al) 



Here M2 is the SU(2) soft breaking gaugino mass, fi is 
the supersymmetric Higgs mixing parameter, tan/3 is the 
ratio of the vacuum expectation values of the two Higgs 
doublets and Mw is the mass of the W boson. 

Similarly the spin- 1/2 superpartners of the W'^ and 
B gauge bosons as well as of the two CP-cven neutral 



Higgs mix after electroweak symmetry breaking. The 
4x4 mixing matrix is given in the bino, wino, Higgsino 
basis by [82| 



/ Ml 






M2 
Mzc^cp 



-Mzs^cp 
Mzc^ cos 13 




-Mzc^sis 




(A2) 

Here Mi denotes the supcrsymmetry breaking bino mass. 
Furthermore = sin 0^ , c„ = cos 0„ and 0w is the 
electroweak mixing angle. Mz denotes the Z boson mass. 
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